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Abstract. We study the spectral properties of a many-body system under a regime 
of electromagnetically induced transparency. A semi-classical model is proposed to 
incorporate the effect of inter-band interactions on an otherwise single-body scheme. 
We use a Hamiltonian with non-Hermitian terms to account for the effect of particle 
decay from excited levels. We explore the system response as a result of varying the 
interaction parameter. Then we focus on the highly interacting case, also known as the 
blockade regime. In this latter case we present a perturbative development that allows 
to get the transmission profile for a wide range of values of the system parameters. 
We observe a reduction of transmission when interaction increases and show how this 
property is linked to the generation of a strongly correlated many-body state. We 
study the characteristics of such a state and explore the mechanisms giving rise to 
various interesting features. 
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The study of light-matter interaction has experienced a pronounced development 
over the last few decades, both experimentally and theoretically. In this process, progress 
has been made by advancing the foundations laid by pioneering observations, such as 
electromagnetically induced transparency (EIT), stimulated Raman adiabatic passage 
(STIRAP) and coherent population trapping (CPT). Such techniques can be understood 
in terms of single particle models, making it possible to obtain much insight into the 
physical background leading to the realisation of the mentioned phenomena. Similarly, 
due to the increasing level of control possible in low-temperature physics, the transition 
toward more collectively driven scenarios is starting to attract general attention. 
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Figure 1. Diagram of EIT including several particles. The dynamics is determined 
by the coupling and probe lasers which drive particles between the ground and excited 
bands and the excited and interaction (Rydberg) bands respectively. In addition, the 
first excitation levels of a system with 8 bosons are shown. 



One way of understanding the mechanism underlying EIT is to think in terms of 
cancellation of excitation pathways with opposite phases [T] . As the process involving 
the absorption of one photon is on-resonance, a different transition enhancing photon 
emission arises with equal intensity and opposite phase. The result can then be seen as 
pure quantum interference. An alternative way is to view the optical field as interacting 
with just one of two possible coupling modes [2]. The resonant mode is related to 
the so-called dark state, which displays a decay time much lower than other possible 
configurations, which allows for the whole system to be pumped into this dark state 
by the combined action of probe and coupling lasers. These views are useful when 
many-body effects are taken into consideration, but signatures of collective behaviour 
slowly degrade the single particle response [3], and a more detailed analysis becomes 
necessary. This response results as a consequence of the interaction among particles and 
is an instance in which photons communicate using matter as an intermediary [1]. In 
this direction, an interesting problem is the study of the many-body physics behind the 
interference mechanism governing EIT. An example of this is the realisation of clean 
EIT-profiles on ensembles of interacting atoms [5l [6] . 

Atoms with very high principal quantum number — also known as Rydberg atoms- 
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display important properties that make them useful in several applications. In 
particular, the interaction among Rydberg atoms gives rise to a blockade mechanism 
originated by dipole-dipole interactions or van der Waals forces. The blockade process 
can be exploited in diverse applications, such as quantum gates proposals [TT | |8| [Tf [9] . 
cavity QED constructs [111 [13], the generation of entanglement [HI [151 IB], adiabatic 
passage [T7|, HE] and numerous quantum information techniques [T9] . 

For the specific case of EIT, inter-atomic interactions produce a highly robust state 
with a collective character. As discussed in reference [5], the quantum superposition 
can be seen as involving states where due to dipole-dipole interactions only a single 
atom can be excited to the interacting band and therefore only one atom is involved 
in EIT. This results in a reduction of transparency as the possibilities of interference 
of excitation pathways are diminished. Such reduction shows neither resonance-shift 
nor line-width broadening. Some of the experimental results cannot be explained via 
mean-field approximations and full quantum models must be employed to reproduce the 
observed features In reference [20] it was shown that two-photon correlations lead 
to an enhanced attenuation of the probe beam for strong intensities and in this way a 
detailed description of the experiment in reference [5] is achieved. When emphasis is 
made on the propagation of the probe laser through the medium as in [21], it is shown 
that the blockade mechanism gives rise to a highly non-local response in addition to 
non-linearities. 

Many-body effects in an ensemble of interacting atoms can be studied using the 
master equation formalism [22]. In this case the number of coefficients necessary to 
describe the density matrix is proportional to the square of the total number of elements 
in the Hilbert space. Instead, we probe the advantages of introducing decay factors 
as imaginary elements in the Hamiltonian, so that the whole analysis can be carried 
in terms of state functions, providing insight into the development of the many-body 
scenario as well as the statistical effects that arise due to the bosonic nature of the 
particles. Non-Hermitian Hamiltonians (NHH) have proved useful in several studies, 
e.g., single particle models of EIT P, STIRAP in the presence of degenerate product 
states [231 l2l] and the dynamics of Bose-Hubbard dimers [25] among many others. This 
approach is reasonable since we want to explore the case of high number of particles 
where standard approximations have displayed mixed results [26]. Below we show 
that the proposed methodology produce consistent results and is especially suitable 
to develop a perturbative approach which is valid over a wide range of parameters. The 
insight acquired in this way is less accessible using a fully numerical approach since 
several tunable parameters must be considered. 

In our proposal we assume that probe and coupling lasers induce particle exchange 
among three energy levels as in a ladder EIT scheme [5l [T7j (figure [1]) and that two- 
body interactions take place only in the interaction band. Such two-body exchange 
is mediated by a constant parameter U . In general, the interaction depends on the 
inter-atomic distance and the Rydberg principal quantum number [5l[T9]. Using a semi- 
classical approach and incorporating the rotating wave approximation in addition to 
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non- Hermit ian decay-terms the NHH reads, 

H = -^]ce~*'^^*a^Jae - fipe~^'^''*ajae + H.c. 

+ ^a^RaRia^RCiR - 1) - ir^ajja/j - iVeCilae 

+ Eghlhg + End^j^ciR + EeCilcie. (1) 

Creation and annihilation operators introduce boson-hke statistics through the 
commutation relations [di, a^] = 1 for L = g,e, R. Qp and Qc are the probe- and 
coupling-Rabi frequencies respectively. Similarly, Up and Uc are the laser frequencies 
while Eg, Er and Ee represent the energies of the corresponding bands. The intensity 
of particle decay from the interaction- and excited-bands is controlled via Tr and Fe 
respectively. Both constants are positive definite. This approach is valid in the limit 
of weak coupling between the ground and excited bands pQ. In our analysis we set 
h = 1 and measure all the parameters in recoil energies Er = h^K'^/{2m). The energy 
structure of the system correspond to three-level EIT-picture where Er > E^ > Eg 
(figure [1]). Although not explicit in ([1]), the total number of particles is M, which 
accounts for the number of atoms in a blockade sphere. The proposed scheme can be 
realised by projecting counter-propagating lasers onto a cloud of ultra-cold atoms — 
for example, ^''^Rb. These lasers provide the probe- and coupling-frequencies in our 
proposal. The number of atoms can be controlled by pumping atoms into energy levels 
that do not couple to the probe- or coupling- lasers. The energy spectrum of the gas 
contains bands that can stand for the ground and excited bands of figure [H The gas 
also contains highly excited atoms which display interaction intensities much larger than 
atoms in the ground or excited levels, so that it is valid to neglect interaction effects on 
these bands. The intensity of the interaction can be tuned using Feshbach resonance 
|27j . Actual experiments implementing this approach have been reported in several 
works, as for instance in references [SI |6] . As a result of the inclusion of non-Hermitian 
terms in equation ([1]), the wave function norm is no longer preserved and therefore the 
quantum state must be normalised in anticipation to any explicit calculation. After an 
appropriate transformation we get the following dressed NHH in the interaction picture, 

Hd = - flc{a''jide + dldn) - Qp^d^de + dldg) 
+ - 1) - i^RaRdji - iTedlde 

+ 5d\dg + {Er - ujc)a^j^dR + Eed^de, (2) 

where 5 = Eg — Up. For simplicity we have assumed single photon resonance and without 
loss of generality we set Er — Uc = Ee = 0. Here we are mainly concerned with the 
magnitude of the coupling between the atoms and the probe laser, i.e., the atomic 
susceptibility. Hence we focus on the mean value, 

Xm - ■ y-^l 

Imaginary and real parts of Xm account for different orders of absorption and 
refraction respectively. Figure H] presents the behaviour of Xm various parameters 
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Figure 2. Xm = 50, = 10, r^, = 0.5, fie = 1, fip = 0.5. Inset: Same 

parameters and 5 — 0. Left: Xm against U . Right:xM against M. The large dip in 
the main figure features a decrease of absorption of the probe laser at zero detuning. 
Increasing the interaction provokes a reduction of transmission. As shown in the 
insets, this effect can be controlled by tuning either the intensity of the interaction or 
the number of particles. 

are tuned. In all cases the state of the system is given by the right eigenstate of equation 
(j2]) corresponding to the eigenvalue with the highest imaginary part (the less decaying 
state). The NHH ([2]) conserves the total number of particles so that the full dimension 
of the Hilbert space is (M + 2)(M + l)/2. The pattern shown in figure [2] is proportional 
to the transmission profile of the probe laser. It indicates a characteristic window of 
transparency that results from the combined action of the incident lasers. In the absence 
of interaction the dip in absorption can be ascribed to interference of pathways with 
opposite phases, i.e., the process by which one atom goes from \g) to |e) is cancelled 
out by the process by which the atom goes from |e) to and then all the way back 
from \R) to |e) to \g) again. As can be seen, the case f/ = 0, which is equivalent to the 
single particle case, displays the maximum interference. As the interaction is gradually 
turned on, the two-photon resonance is shifted to the right due to the energy increase of 
the interaction band produced by the repulsion among particles, xm continues to grow 
until it reaches a saturation value, but without completely suppressing transparency. 
Xm also grows with M, but no saturation is visible over values of M less than 100. 
In general, we can see that the insertion of particles causes a reduction of quantum 
interference. Such an effect can be enhanced either by increasing U or M. However, in 
every case the physical response is different. While changing U affects the intensity of 
two-body interactions, adding particles produces a (sometimes steep) rearrangement of 
the quantum state. 
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From the Heisenberg equations we can extract the expression, 

da\ t + 

iVLpOil - i5dr (4) 



where a\ = e~**-'^-DQ,t git-f^u ^qj- ^ — g^^^ g_ gjnce VLp « Qc we can assume that does 
not greatly influence the evolution of and dj, so that al can be treated as a function 
of time. Thus it follows that, 



''5 



&l{t) = iQpC-'^' / e'^^al{T)dT. (5) 



Under such an assumption, in the blockade regime equation ([21) can be reduced to 
a Jaynes-Cummings-like Hamiltonian on the interaction- and excited-bands. 

Hue = -^c{cr^a^ + 0-^4) - "^^R ~ 

where we have employed Pauli matrices to account for the dynamics of the 
interaction band. The eigensystem of equation ([6]) is given by [28], 



E = 0, |0,0), (0,0|, (7) 

E^ = -t{T,n + Tn + n^V^T+Te^''") , (8) 

\E^) = {\l,n)+ze^''-\0,n+l))/Zt (9) 

{E^\ = {{l,n\+ze^''-{0,n + l\)/Zt (10) 

where, 

z± = yir;^, (11) 

cosg„= n = 0,l,...,M-l. (12) 

The first and second integers in a ket (bra) correspond to the number of particles in 
the interaction (excited) band. The angle 6n is complex as well as the eigenvalues of the 
NHH. In the form in which they appear above, the eigenvectors satisfy {El\E^,) = S^'S^' . 
This eigensystem can be used to write, 

M-l 

Hr, = Y.E^\Ei){E^\. (13) 

n=0 

The spectrum of the NHH is well defined except when cos On = 1, in which case the 
eigenvectors ^ and flTUl) become indeterminate but one can still recover the NHH as the 
limit of equation (fT3|) . These singularities therefore produce removable discontinuities 
leading to no divergence in the mean values of the system's observables. Now we project 
on the subspace associated to the NHH, 

at = at r^Vn)(^n I + |0,0)(0,0|| . (14) 

I n=0 ) 

This procedure is facilitated by introducing the coefficients validating the following 
identities: 
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al%Q) = ct\E^) + c,\E,) , (15) 



«Il^n ) = Ca+l\En+l) + C„+i|^„,+i), (16) 
l\En) = dn+l\Ki+l) + dn+l\En+i). (17) 



Once a], is written in terms of the eigenvectors of the NHH we can calculate a\ and 
replace it in equation ([5]) where we can now carry on the integration in terms of time. 
The result can be written in the form, 



a\{t) = al + np(k{f))-e-'''k{t)), [li 



where, 



+ E^5 + 

Let us assume that time- dependent terms decay over time. The range of validity of 
this assumption will be addressed further below. Since at t = all the particles remain 
in the ground band the stationary state can be obtained from in accordance to the 
Heisenberg picture, 

{alioo)f 10) = E ( ^ ) {^if'" i^pmriO)- (20) 

Operator k{0) couples the representations of Hrs in such a way that, 

k{0)'\0) = v,ZU\Et,) + w,Z,_,\E^_,),k = 1, . . . , M. (21) 
By replacing in (I2U]) one obtains. 



M 



M 



+i{vke''''-' + Wke-'''^~')\M - k, 0, k)) . (22) 
The set of coefficients Vk,Wk are connected by the recursion expression, 

f "0=^^f 1' k = 2,...,M. (23) 
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Figure 3. First order absorption (a) and second order refraction (b) profile for different 
M . The cases M — 1 and M — oo correspond to equations (1551) and ([59)1 respectively. 
Other cases shown are generated following the procedure sketched in the text. The 
parameters involved in the computations are Fg = 2, F^j — Q, VLp = 0.1 and fie = 1. 



From an induction argument is found to satisfy, 
Bk2i\/ksin6k-i = 



(24) 



-ie 



(i?^_i-Sfc"_2-5) 



■5) 



■5) 



and the first elements are given by, 

_ -1 _ 1 

~ 2(^0+ -5) sin 00 ' ~ 2(^0- -5) sin^o' 
Equations (!25|) . (El]), (!23|) and (!22|) can be used to numerically generate the 
stationary state for any set of parameters excluding singularities. Such a condition 
can be improved by introducing new discrete variables, 

Pk = Vk + Wk, qk = sm6k-i{vk - Wk). (26) 

These variables obey a recursion relation analogous to equation ( ]23|l with Vk and 
Wk replaced by pk and qk respectively and Bk replaced by. 



A, 



in such a way that. 



1 / RkUJ + SkTk sin^ 9k^2 RkTk + Sk(^ 



y/kFk \ rkUJ + SkTk sin^ 9k-2 TkTk + SkUJ 



Fk = uj^ - Tl sin^ dk-2, 



(27) 



Many-Body Effects in a Model of Electromagnetically Induced Transparency 



9 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 




1 

0.5 




Regular Relaxation 



Exponential 
Growth 



1000 



2000 3000 4000 5000 6000 7000 8000 9000 10000 
M 



Figure 4. 
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In every case (5 = 0, F/f = Q, VLc = 1 and ilp = 0.1. The 



absorption curves display non-linear signatures. Saturation takes place according to 
the interplay between decay and coupling. Inset: Relaxation map for M > \. The 
case M — 1 shows no signatures of exponential growth. 



Rk = i^cJ cos ^fe-i - Wk - 1, (2^ 



Sk = -n,{2k-l), (29) 
rk = _ n,^k{k-l) (sin^ + sin^ Ok^^) , (30) 

Sk = i^c cos 9k-i ~ rjVk, (31) 
Tk = ~2r]n,Vk - 1, (32) 

77 = 2!^ + S. (34) 
These expressions along with the initial coefficients, 

(zr^ + 5)(^re + 5)-fir i-^-^^ 

can be integrated into a programming routine that recursively calculates the state 
coefficients and then xtd ■ We note that when — TrI << 2Qck the recursion matrix 
can be approximated as follows, 
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Employing this one finds a stationary state displaying the following susceptibility 



Xoo / '"n I r\ ' Xoo (Xoo ) * 



Similarly, from a direct calculation the single particle susceptibility is found to be, 

(1) ^ -%{l^R + 5) , . 

(,r^ + 5)(,r, + 5)-fir ^ ^ 

As can be seen in figure [H equations ( !38l) and ( !39l) are extreme cases of the numerical 
results obtained using the matrix recursively. Here we have chosen the case = 
because it is most close to the actual experimental case where the decay rate of the 
Rydberg state is usually negligible. The inclusion of more particles in the system induces 
a rise of Im (^Xm ) ^^ro detuning. In the same instance. Re (xi? ) undergoes a change 
of sign as it dips from its peaking value. From figure H] we can see that the absorption 
profile at 5 = features non-linear behaviour over a wide range of values of M and then 
asymptotically converges toward the estimation given by equation (138|) . As we increase 
Fe, the growth speed decreases and the curves display stronger inflection. In every case 
the absorption profile is characterised by a steep absorption growth. Such a non-linear 
response is linked to a cooperative many-body state where the single particle outcome 
is no longer dominant. As M goes up further, the absorption value starts to saturate, 
suggesting that less particles are being integrated in the interaction process. 

In the single particle case, maximal transmission is achieved as a result of almost 
perfect quantum interference between excitation pathways. In this case the only coupling 
mode involved in EIT in equations (l8l fT2|) is n = 0. Gradual increase of Vlp allows more 
particles into the excited band and therefore more coupling modes participate in the 
EIT process. One important characteristic of EIT is that the reduction of absorption 
at zero detuning is accompanied by a peak in the refraction coefficient [1] . This allows 
for powerful applications such as slow light or light storage. The mechanism behind 
the latter procedures is based on the fact that at 5 = absorption can be made low at 
the same time that the coefficient that determines the velocity of light in the medium 
peaks. This important property can be observed without much effort using our approach. 
Figure [3^ shows a dip in Im (^Xm ) while figure [3]d shows and a peak in Re (xm ) the 
same parameters. Also at single-photon resonance it can be shown that independently 
of M, 

^''i?^M),^, = I^(x^^}),^, = ^, (40) 

(2) 

Notice that the peaking of refraction can only be seen for M > 1 because Xi 0- 
This means that the system's refraction at 5 = depends only on two-or-more-body 
transitions, where several atoms are excited (or decay) simultaneously. In a sense, one 
can think of the process taking place at zero detuning as one in which single-body 
transitions play a rather marginal role and instead the leading response is mediated 
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by higher order transitions. Likewise, couphng modes involving more than one particle 
display less interference, as many particle pathways between the ground- and excited- 
bands can only interfere partially with their blockade counterparts. 

As even more atoms are integrated in the model, the possibilities of light being 
absorbed via particle excitation become higher and the characteristic profile of EIT 
finally fades away. The enhancement of absorption due to many-body effects is especially 
notorious in the range M < 1000. Since the peak of refraction at M = 2 is positive one 
can find a value of M (M = 65 for the values of figure [3]) for which refraction is almost 
zero at zero-detuning, while absorption is still low. In this case the system is almost 
unresponsive to the incoming radiation. 

In order to check for the consistency of our approach we stress that our central 
assumption is the cancellation of k{t) in equation ( fT9l) as t — )■ oo. This is indeed the 
long as all the imaginary parts of the arguments of the exponentials in equation 
(fT9|) turn out to be negative. Only the following arguments could become greater than 
zero, 

/m(i?„V - E-) = 1 {D{n) + D{n + 1)) - T,, (41) 

and, 

ImiE-^, - E-) = 1 (Din) - D{n + 1)) - T,, (42) 

where. 

Din) = Re i^^ [T, - Tr)^ - AQl{n + , (43) 
for n = 2, 3, ■ ■ ■ , oo. 

If either (14T|) or (H2|) become positive then operator k{t) will display exponential 
growth and will dominate the stationary state. While it seems operationally possible to 
obtain such a state, this feature is less consistent with our initial consideration in which 
Qp is a perturbative parameter and hence most particles remain in the ground band. 
Due to the form of D{n), the maximum value of ( HTl) and ( H2l) takes place at n = 2. 
Therefore, if either (HT]) or (H2|) are positive for any n > 2 then they are positive for 
n = 2 as well. Hence n = 2 is the only relevant mode in a relaxation analysis. We have 
depicted in the inset of figure H] the relaxation map that results following the arguments 
discussed above. It is worth pointing out that any arbitrary set of realistic parameters 
in which > fall well inside the regular relaxation zone. It also becomes apparent 
that in the many-body case the mere existence of a dark state does not guarantee the 
state convergence toward such a state. 

We have investigated the reduction of transparency in an EIT set-up as a result 
of the collective character developed by many-body matter. Results corresponding to 
a semi-classical model were obtained from numerical diagonalisation as well as from a 
perturbative approach. In the latter case we presented a semi-analytical procedure that 
can be used to find the stationary state of the system in the limit of small intensity of the 
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probe laser. We have also studied the range of validity of our method and established 
explicit conditions for regular relaxation. The procedure itself shows interesting issues 
and is valid in the range of realistic experimental parameters. We have in this way 
presented an alternative analysis of the effects of many-body interaction on EIT. As a 
prospect extension of the present work, it would be interesting to introduce a light mode 
to describe the probe laser in order to explore the evolution of initially-coherent states 
of light and the effect of particle interaction on photons. 
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